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' The Casimir problem is usually posed as the response of a fluctuating quantum field to externally 

' imposed boundary conditions. In reality, however, no interaction is strong enough to enforce a 
. boundary condition on all frequencies of a fluctuating field. We construct a more physical model of 
I the situation by coupling the fluctuating field to a smooth background potential that implements the 
I I boundary condition in a certain limit. To study this problem, we develop general new methods to 
compute renormalized one-loop quantum energies and energy densities. We use analytic properties 
• of scattering data to compute Green's functions in time-independent background fields at imaginary 
^ . momenta. Our calculational method is particularly useful for numerical studies of singular limits 
^ , because it avoids terms that oscillate or require cancellation of exponentially growing and decaying 
H ' factors. To renormalize, we identify potentially divergent contributions to the Casimir energy with 
- - ' low orders in the Born series to the Green's function. We subtract these contributions and add 
back the corresponding Feynman diagrams, which we combine with counterterms fixed by imposing 
standard renormalization conditions on low-order Green's functions. The resulting Casimir energy 
and energy density are finite functionals for smooth background potentials. In general, however, 
the Casimir energy diverges in the boundary condition limit. This divergence is real and reflects 
the infinite energy needed to constrain a fluctuating field on all energy scales; renormalizable 
quantum field theories have no place for ad hoc surface counterterms. We apply our methods to 
simple examples to illustrate cases where these subtleties invalidate the conclusions of the boundary 
condition approach. 
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1. Introduction 



Long ago, Casimir predicted the existence of a force generated by quantum fluctuations of 
the electromagnetic field between uncharged, perfectly conducting metal plates in vacuo |1| . 
Since then, many related calculations have been carried out (for reviews see ||, |, ||). 
The standard approach to problems of this type has been to consider the response of a free 
quantum field to time-independent boundary conditions imposed on stationary surfaces. 
This calculation is an idealization of the more physical description of the quantum field as 
interacting with a smooth external potential, which goes to zero away from the surfaces. 
This paper will provide the technical tools necessary for understanding precisely when this 
idealization is justified, and when it is hazardous. We will demonstrate the use of these tools 
in some simple models; the reader is also referred to [0] for a briefer, less technical, and more 
focused discussion of a variety of Casimir problems. 

In the field theory picture, the dynamics of a fluctuating field in a background potential 
is governed by the one-loop effective energy, which can be expressed either as an infinite 
sum of Feynman diagrams or as a "Casimir sum" over modes. Both expressions are formally 
divergent, so they can only be defined by introducing counterterms fixed with physical renor- 
malization conditions. A Dirichlet boundary condition emerges when a repulsive background 
scalar field becomes strong and sharply peaked. We will see that the renormalized total en- 
ergy of the fluctuating field relative to its value in the absence of boundaries can diverge in 
this limit. In such cases, the idealized Casimir energy per se does not exist. However, these 
divergences are associated with the boundaries and therefore do not affect the energy density 
away from the boundaries or the force between rigid surfaces. Nevertheless some geometric 
quantities, such as the "Casimir stress" |Q — the force per unit area on the surface where 
the boundary conditions are applied — are sensitive to these divergences. In such cases, 
there is no way to calculate meaningful finite results by considering boundary conditions. 

To demonstrate these results in detail, we require efficient and robust techniques for 
computing Casimir energies and energy densities, which we describe in Section 2. The 
calculation of Casimir energy densities is interesting in its own right, for example in the 
study of energy conditions in general relativity |^, |TU]. Our approach will phrase the 
calculation in the language of conventional renormalizable quantum field theory, in a way 
that is amenable to efficient numerical computation. We study the vacuum polarization 
energy and energy density of a fluctuating boson field coupled to a background cr{x), 
which is sharply peaked in the boundary condition limit. We introduce counterterms fixed 
by perturbative renormalization conditions, which define physical inputs to the theory such 
as particle masses and coupling constants at particular values of the external momenta. The 
resulting renormalization scheme is conventional, precise, and unambiguous. It can be related 
to any another conventional scheme by the analysis of low-order Feynman diagrams. We do 
not, however, see any relation between this scheme and analyses in which boundary surface 
dependent counterterms are introduced in an ad hoc manner. It is a common procedure in 
the literature to simply ignore divergent surface integrals and pole terms arising in the course 
of the calculation |Tl|, O]. The (^-function method, in particular, introduces counterterms of 



any power of the geometrical parameter |]13l . Such counterterms redefine classical parameters 
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of the background, but have no justification in quantum field theory. As a result, these 
approaches do not appear to be reconcilable with ordinary renormalizable quantum field 
theory. 

Our method is based on an extension of methods developed for the total energy [|14|, [15[| , 



together with an extension of the analytic continuation methods introduced in Ref. p6[| . 
First, we relate the matrix element of the energy density operator to the scattering Green's 
function at coincident points, G{x, x, k). We regulate the large k behavior of G by subtracting 
the first few terms in its Born expansion and adding back the contribution of the associated 
low-order Feynman diagrams [|1^, |T^. The Green's function is difficult to deal with at large 
(real) k when the background field approaches the boundary condition limit: it oscillates 
with amplitude much larger than the net contribution to the energy density. To avoid such 



oscillations we rotate k to the imaginary axis |16, 17|. Before doing so, however, we rewrite 



the Green's function as a product of terms, each of which is bounded for k on the imaginary 
axis and each of which can be easily computed by integrating a Schrodinger-like equation 
subject to simple boundary conditions. This formalism allows us to avoid exponentially 
growing quantities anywhere in the calculation. As a result, we are able to write the Casimir 
energy density as an integral along the positive imaginary axis and perform the integral 
directly for imaginary k. Here we depart from Refs. |T^, |T^, which formulate the problem 
as an integral along the imaginary axis but then use a dispersion relation to compute the 
integrand in terms of scattering data for real k. 

The final ingredient we need for the calculation is the Feynman diagram contribution, 
which is where we implement the renormalization conditions. The energy density operator 
Too generates insertions in Feynman diagrams of every order in the background field a. Using 
a functional approach allows us to keep careful track of orders in a in this calculation, which 
we need to maintain consistency with the rest of the calculation. 

In Section 3 we show how to apply our methods to a simple example, the Casimir energy 
density for a boson field in one space dimension. We couple the fiuctuating field to a 
background (y{x) with an interaction of the form A</)^ct(x). For (j{x) = 6{x — a) + 6{x + a), we 
obtain Dirichlet boundary conditions at a; = ±a in the limit where the coupling constant A 
becomes infinite. We find that the total Casimir energy diverges in this limit. The divergence 
is localized at the points where the background a is nonzero. For all x ^ ±a, the energy 
density remains finite and as A approaches infinity it smoothly approaches the result obtained 
by imposing the boundary condition = at x = ±a. To complete Section 3 we study the 
vacuum polarization energy in the presence of a smooth, strongly peaked background field. 
By taking a to be the sum of GauBians centered at ±a, we can study the way that the 
energy density behaves as a{x) approaches the boundary condition limit. As expected, it 
approaches its limiting form at every x ^ ±a, though the rate of approach depends on the 
proximity to ±a. 

In Section 4 we consider the case of a circle in two spatial dimensions. We consider 
GauBian backgrounds peaked around r = a and study the limit in which the GauBian 
approaches a 5-function. When the strength of the coupling goes to infinity, we obtain a 
Dirichlet boundary condition at r = a. However, the renormalized total energy diverges 
in the 5-function limit, even at finite coupling. This result is a simple consequence of the 
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fact that / (Px[a{r)Y' — > oo in this hmit. As a result, the surface tension on the circle — 
the derivative of the energy with respect to a — diverges. This example illustrates the 
generic problem we have found with attempts to define the Casimir energy in the context of 
a quantum field theory: Even though the renormalized vacuum polarization energy is a finite 
functional of a smooth background, it diverges as the background field assumes the singular 
configuration necessary to implement the boundary condition. This divergence affects any 
physical quantity, such as the surface tension, whose measurement requires comparison of 
configurations with different boundaries, for which the divergent contributions are different. 
We take this divergence as an indication that such quantities depend on the details of the 
interaction between the fiuctuating field and the material to which it couples. Thus the 
Casimir stress cannot be defined in a way that is independent of the other stresses to which 
the material is subject. 

We close with a summary and some considerations for future work in Section 5 and give 
technical details in the Appendices. 

2. Method 



We begin by developing the method for computing the energy density of a fiuctuating 
quantum field coupled to a classical background. Prom the corresponding scattering and 
bound state wavef unctions, we construct a Fock decomposition of the quantum field and 
compute the matrix element of the energy density operator Too- We express this matrix 
element in terms of a Green's function with appropriate boundary conditions. The energy 
density is given by a sum over bound states plus an integral over the continuum labeled by 
the wave number k. We develop a representation of the Green's function suitable for analytic 
continuation into the upper half /c-plane, so the k integral can be deformed along a cut on 
the positive imaginary axis. To regulate the ultraviolet divergences of the theory, which cor- 
responds to eliminating the contribution associated with the semi-circle at infinite complex 
momenta, we subtract leading Born approximations to the Green's function, which we later 
add back in as Feynman diagrams. These diagrams are then regularized and renormalized 
in ordinary Feynman perturbation theory. 

2.1 Formalism 



We consider a static, spherically symmetric background potential a = a(r) with r = l^l 
in n spatial dimensions. The symmetric energy density operator for a real scalar field coupled 
to a is 



Too{x) 



0' + (V0)' + m^c/)^ + (7(r)02 
0' + 0(-V' ' 



m 



a(r 



where we have rearranged the spatial derivative term in order to be able to use the Schro- 
dinger equation to evaluate the expression in brackets. We assume that the background 
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potential is smooth, with dr\a{r)\{l + r) < oo.0 We define the "vacuum" to be the state 
\fl) of lowest energy in the background a and the "trivial vacuum" to be the state |0) of 
lowest energy when a = 0. The vacuum energy density is the renormalized expectation value 
of Too with respect to the vacuum (fi|Too(x)|f2)ren, which includes the matrix elements 
of the counterterms. Since we have spherical symmetry, we will define the energy density in 
a spherical shell, 

<r) = j^r--'{moix)\n),,n. (2) 
We first perform a partial wave decomposition using spherical symmetry, 

<f){t,x)=Y.Mt,r)y{ey{x) (3) 

where y{e}{x) are the n-dimensional spherical harmonics{£} refers to the set of all angular 
quantum numbers in n dimensions. The total angular momentum assumes integer values 
£ = 0, 1, 2, ... in all dimensions except for n = 1, where i = and 1 only, corresponding to 
the symmetric and antisymmetric channels respectively. We make the Fock decomposition 

Mt,r) = ^ / -y= \Mk,r)e-''^'a,{k)+^;{k,r)e''^'al{k) 

r~2~ Jo y/7^UJ L 

+i E [M^) + Mr) , (4) 



f 2 
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which we have split into scattering states with u = \^W+m? and bound states with uj 



m? — nj. Note that the latter wave functions are real. The vacuum is annihilated by all 
of the ai{k) and a^j. The radial wavefunctions in eq. (§) are solutions to the Schrodinger-like 
equation 

-< + l(z/-i) (z/ + i) ^ + a(r)^-A;2^ = (5) 

where 



Tl 

u = i-l + -. (6) 
In each channel, we normalize the wavefunctions to give the completeness relation 

2 roo 

- dk^*^{k,r)^g{k,r') +y^4jij{r)^(>j(r') = 5{r -r') (7) 
which implies the orthonormality relations 

roo 

/ ipij{r)^piji{r) dr = 6jji for the bound states and 

/•oo ° 

/ tpUk,r)tjj£{k' ,r) dr = —5{k — k') for the continuum states. 

Jo 2 



^In view of later applications with singular background fields, we can relax this condition to allow for 
all a{r) for which a scattering problem in the usual sense can be defined. However, some of the standard 
bounds on the asymptotics of Jost functions no longer hold for a non-smooth (T(r), giving rise to the 
additional singularities discussed later. 
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For cr = 0, the normalized wavefunctions are il)f'\k^r) = ^J^krJi,{kr). With the above con- 
ventions and normahzations, the standard equal time commutation relations for the quantum 
field yield canonical commutation relations for the creation and annihilation operators, 
[a(,{k),a\,{k')] = 5{k — k')Si/' and [a£j,al,j,] = 5jji5u> while all other commutators vanish. 

The vacuum expectation value of the energy density can now be computed by using 
eq. in eq. (|I]). We obtain 



dk 



—ujilj*^{k, r)^e{k, r) + ^ '^^j. 



IT 



Jo nu ■ zUi 



e(o)(r) + ecT(r) 



(9) 



where Nf 



[n 



2i — 2) is the degeneracy factor in the i partial wave, D 



r{n-i)T{e+i) ' 

'§r ~ ^^)' ^CT{f) is the counterterm contribution, and e*^°^(r) indicates the subtraction 
of the energy in the trivial vacuum, which is just given by evaluating the first integral 
with free wavefunctions il)f\k,r). This subtraction corresponds to a renormalization of the 
cosmological constant. 

We can identify the scattering state contribution with the Green's function by defining 
the local spectral density 

k 

pe{k,r) = -Gi{r,r, k) 



in the upper half-plane, so that for real k 

Re{pi{k,r)} = ijj}{k,r)tjje{k,r) = \m{kGi{r,r, k)} 

where 



(10) 



(11) 



Ge{r,r',k) 



TX JO 



- / dq 



{k + leY — 



E 



k^ + K^, 



(12) 

The ie prescription has been chosen so that this Green's function is meromorphic in the 
upper half-plane, with simple poles at the imaginary momenta k = iKj corresponding to 
bound states. Thus we have 



°o dk 



2m 



-u 



1 + 



4cj2 



1 + 



Dr kG£{r,r,k) 

7A,,-(r)2-e(0)(r) + ecT(r) 



(13) 



where we have made use of the fact that for real k, the imaginary part of the Green's 
function at r = r' is an odd function of k, while the real part is even. We would like to 
use eq. (p!3|) to compute e(r) as a contour integral in the upper half-plane. First we must 
eliminate the contribution from the semi-circular contour at large \k\ with Im(fc) > 0. Using 
the techniques of Refs. |18| and [O] we see that subtracting sufficiently many terms in the 
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Imk 



iK|i 

iK2, 

iK3i 



Rek 



Figure 1: Contour for k integration. 



Born series from the Green's function yields a convergent integral. We then add back exactly 
what we subtracted in the form of Feynman diagrams. We define 



[piik,r)]N = Pe{k,r) - pf\k,r) - pY\k,r)... 



{k,r] 



k r 



Ge{r, r, k) - (r, r, k) - G\'\r, r, A;) - ... - Gy'^(r, r, A;) 



(14) 



where the superscript (j) indicates the term of order j in the Born expansion. Subtracting 
the free Green's function Gf\r, r, k) corresponds to subtracting e^^\r) above. We substitute 
[pi{k,ry\j^ for kG£{r,r, k)/i in eq. (0) and add back in the Feynman diagrams correspond- 
ing to the subtractions in eq. (0), which we compute below. As a result, we obtain a finite 
expression amenable to contour integration. Furthermore, we have precisely identified the 
potentially divergent pieces with Feynman diagrams, which we can regularize and renormal- 
ize using standard methods. When combined with the contribution from the counterterms 
^CT^r) they yield finite contributions to the energy density (for smooth backgrounds). 
We now have to specify the contour along which to compute the integral 



Jc 27r 



[peik,r)] 



N 



(15) 



The integrand has a branch cut along the imaginary axis k G [im, +ioo] and simple poles at 
the bound state momenta k = iKj. We deform the integral around the cut as shown in Fig. 0, 
picking up the residues, which cancel the bound state contributions in eq. (|13D |Q. Note 
that the Born terms do not introduce any poles. Hence we are left with the discontinuity 
along the cut. 



e r 



°o dt 



TT 



4(t2 -m2)^' 



N 



i=l 



-fdI 
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N 

^Kr) + T.4l{r) + ecT{r). (16) 
1=1 

Here epj)(r) is the Feynman diagram contribution at order i in the background field. As 
emphasized in the Introduction, the representation of e(r) as an integral along the imaginary 
axis is only useful if it can be computed efficiently. In the following subsections we show 
how to do so. 

2.2 Computational Techniques 



In this subsection we concentrate on e(r), the continuum integral in eq. (jlGl) . We develop 
a representation of the partial wave Green's function, Gg{r, r', k) and its Born series, eq. (H), 
that can be computed efficiently for imaginary k. By working with pure imaginary momenta 
we avoid functions that would oscillate at large \k\. We will make extensive use of the 
analytic properties of scattering data 

We begin by introducing solutions to eq. (^ that obey a variety of different boundary 
conditions. First we define the free Jost solution 



we{kr) = {-ly^l^ [J,{kr) + tY,{kr)] , (17) 

where z/ = £ — 1 + |. The radial function we{kr) is a solution to eq. (|^) with cr = describing 
an outgoing spherical wave. Then we define 

• The Jost solution, fi{k,r). It behaves like an outgoing wave at r — oo, with 

^"3^ — 7rT = ^- 
— Wi[kr) 



r— >oo 



The Jost function, Fe{k). It is obtained as the ratio of the interacting and free Jost 
solutions at r = 0, 

F,(A:) = lim/%4- (^9) 

The regular solution, (f)i{k,r). It is defined by a /c-independent boundary condition at 
the origin, 

r(i/ + 1) /r\-('^+i) 



"iS^Vj^Uj ^'ft--)-!- (2°) 

The physical scattering solution, ipi{k,r). It is also regular at the origin (and thus 
proportional to (f)i{k,r)), but differently normalized^ 

k'^+i 

Mk,r) = jr^^Mk,r). (21) 



■^Whenever a fractional power of k appears, we define it to be the limit as k approaches the real axis from 
above. 
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The reason for distinguishing two regular solutions is that (pi has a simple boundary 
condition at r = 0, while ipi has a physical boundary condition at r — oo, corresponding to 
incoming and outgoing spherical waves. With these definitions, the Green's function has a 
simple representation 

G,(r.r'.t) = M^i^|^lIi)(_,)^4, (22) 

where r> (r<) denotes the larger (smaller) of the two arguments r and r' . The poles of G 
occur at the zeros of the Jost function, which are the imaginary bound state momenta. Note 
that these are the only poles of eq. (0) in the upper half-plane, and since the two functions 
in eq. (^2l) and eq. (p2|) obey the same inhomogeneous differential equation they are indeed 
identical. 

The representation of Gi{r, r, k) in eq. (^) is not yet suitable for numerical computation 
on the imaginary axis. Although is analytic in the upper half-plane, and 0^ contain 
pieces that oscillate for real k and exponentially increase or decrease respectively when k 
has an imaginary part. We are actually only interested in the case r = r', in which case the 
product fi(f)i is well-behaved. But because of the exponential behavior of the two factors, 
the product cannot be computed accurately by simply multiplying the individual functions 
obtained from numerical integration of the respective differential equations. To proceed, we 
factor out the dangerous exponential components with the following ansatz^ 

feik,r) = weikr)gi{k,r) 

and 

Mk,r) = — -, (23) 

lu we[kr) 

where Wi is the free Jost solution introduced above. With these definitions. 

As we shall show, both g^lk, r) and h^lk, r) are well-behaved on the imaginary axis. Note that 
the definition of hi does not just remove the free part.Q Instead, it enforces the cancellation of 
Wi in the Green's function. The factors in eq. (p3|) were chosen to provide a simple boundary 
condition on he{k,r) at r = 0. After analytically continuing to k = it, gi{it,r) obeys 

g'litt, r) = 2tUtr)g'i{it, r) + a{r)gi{it, r) (25) 

with the boundary conditions 

\mi gi{it,r) = 1 and lim (^^(it, r) = , (26) 

''^For n — 1 and n — 2, the case of £ = is somewhat different. We will explore the former in Section 3 
and the latter in Section 4. 

^Removing the free part would correspond to an ansatz like ~ wihi as in computations of functional 
determinants with Euclidean Green's functions [E2|. 
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where a prime indicates a derivative with respect to the radial coordinate r. Using these 
boundary conditions, we integrate this differential equation numerically for gi{it,r), starting 
at r = cx) and proceeding to r = 0. Similarly, hi{it,r) obeys 



h'l{it,r) = -2t^e{tr)h'^{it,r) + 



air 



2r 



dUr) 



dr 



T=tr. 



hi{it,r) 



with the boundary conditions 

he{it,0) = and h'^{it,0) = 1 
and we integrate numerically from r = to r = oo. For real r. 



-^In [weiir)] 
dr 



(27) 



(2J 



(29) 



is real with lim.r-,00 ^e{'^) = 1, so the two functions he{it,r) and ge{it,r) are manifestly real. 
They are also holomorphic in the upper fc-plane and, most importantly, they are bounded 
according to 



\g£{k, r)\ < const 
\hi{k, r)\ < const - 



2z/r 



\k\r 



(30) 



so that neither gi nor he grows exponentially during the numerical integration. Thus the 
representation of the partial wave Green's function in terms of ge and hi is smooth and 
numerically tractable on the positive imaginary axis. We have avoided oscillating functions or 
exponentially growing functions that eventually would cancel against exponentially decaying 
functions. 

Finally, the computation of the Born series, eq. (p!^, is also straightforward in this 
formalism. We expand the solutions to the differential equations eq. (|25| ) and eq. (127|) about 
the free solutions, 



ge{tt, r) = 1 + g^"^ {it, r) + gf> {it, r) + . 
ht,{it,r) = 2vrI^{tr)K^{tr) + hf\it,r) + h\"\it,r) + . . . , 



(2), 



,(2), 



(31) 



where the superscript labels the order of the background potential a. The higher order 
components obey inhomogeneous linear differential equations with the boundary conditions 



lim g^p{it, r) 








and lim g^p'{it,r) - 
and /i?'^'(it,0) = 0. 



0, 



(32) 



In these equations a is the source term for g^^\ ag^^^ is the source term for g^'^\ and so on. 
We then obtain the Born series for the local spectral density by substituting these solutions 
in the expansion of eq. (E3) with respect to the order of the background potential. Thus 



10 



we have developed a computationally robust representation for the Born subtracted local 
spectral density, 

h(,{it,r)gi{it,r) 



N 



2ugi{tt, 0) 



(33) 



N 



2.3 Feynman Diagram Contribution 

In this subsection we develop methods to compute the contribution from the Feynman 



N 



diagrams and counterterms, J2 CfdI'") +^CT{'f'), in eq. (ffO). To one-loop order, the Feynman 

i=0 

diagrams of interest are generated by the expansion of 



|Too(a;)|0)~-Tr 



-1 



(34) 



to order N in the background a. Here is the coordinate space operator corresponding to 
the insertion of the energy density defined by eq. (|l|) at the spacetime point x, and the trace 
includes space-time integration. Since we are doing ordinary perturbation theory, we take 
the matrix elements with respect to the trivial vacuum, which is annihilated by the plane 
wave annihilation operators. The energy density operator has pieces of order a" and cx^. 



00 



d'ym (Tf + Ti^)) 0(2/) (35) 



where d = n + 1 is the number of space-time dimensions. We have separated the different 
orders in the external background field. 



dt S{x — y) dt + V 6{x — y) V +m'^5{x — y) 



t« = a(x)5(x-y). (36) 
The 5-functions are understood as rf-dimensional. The arrows denote the direction in which 



the derivatives act when substituting this representation into the functional trace in eq. ( p^ ) 
and y is the coordinate space position operator, with y\^) = ^|^). However, it is more 
convenient to perform the computation in momentum space. The relevant matrix elements 
are 



(fc'iffiA:) 
(A;'|f«|fc) 



^i{k' —k)x 



k^'k^ + k' -k + m' 



a[x]e 



i{k'—k)x 



(37) 



Here we will explicitly consider the contributions to (Too) that are linear in a because these 
are the only contributions required in the examples studied in Sections 3 and 4. The first 
piece comes directly from Tj^^ as shown in the left panel of Fig. 0, 



-Tr 



1 



— (9^ — 



fa) 



-a(x] 



d'^k 



(3J 
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J- Y 




Figure 2: First order Feynman diagrams contributing to the energy density (Too(x)). The soHd 



circle represents the cf) loop and the operator insertions from eq. (36) are denoted by T^'^\ In 
the right panel, a is the local insertion of the background field from expanding the propagator, 
cf. eq. 



This local contribution is ultraviolet divergent for d > 2. It is canceled identically by the 
counterterm in the no-tadpole renormalization scheme, since the right-hand side of eq. ( p8| ) 
equals minus the counterterm Lagrangian evaluated for the background field a. For further 



discussion of this renormalization condition and its consequences, see Ref . [jT5| . The second 
contribution at order a originates from T^^^ and the first-order expansion of the propagator 
as shown in Fig. 0, 



iTr 
2 



—(9^ — m? 



-d"^ — m? 



(39) 



For this diagram we find 

1 r d'^k j 

2 J {27iY J {27iYk^-m^' 



d^k' 



(k-k' + 



m 



1 



d'^q ^ 



{2n] 



^qx I 




d'^l l-l + m'^ 



_ ^2 

C(i-C)g-g 



a{k - k') 



(40) 



where a{q) = 27T6{q^)a{q) is the Fourier transform of the (time independent) background 
field. We have defined p = {p^, —p) for the Lorentz-vector p = {p^,p). The leading ultraviolet 
divergences cancel in eq. (^) between the temporal and spatial components. We thus obtain 
the expression for the Feynman diagram and counterterm contribution through first order 



m a 



-fdI 



2iT 2 r 



ecT(rJ 



r(f) 



2 



1 



1 



m 



2^x 



— m? 



a 



2t[ 2 r 



I .-ir 2- 



d 



r(f) 



(4vr)d/2 



-^diq)e"^-^ C dC 
{2'kY Jo 



C(i-C)g 



m2 + C(l-C)g' 



2-d/2 



(41) 



We note that this piece does not contribute to the total energy because it vanishes when 
integrated over space. The extension to higher order Feynman diagrams in a is straightfor- 
ward. For example, the contribution O^a"^) is obtained from the two terms 

1 1 \ i„ 1 1 1 



iTr f« a 

2 V ^' 92 + m2 d^ + m^ 

We use the matrix elements, eq. 



- -Tr ( f 



ra- 



rer- 



(42) 



a(g)e*'?^ / dC 



to write this in terms of ordinary loop integrals 

1 



(27r 



{27rY [/2 - m2 + C(l - C)g2]^ 
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i f d'^q f d'^p ^ 



X C dc dn f '^''^ ^■i + [Cq + vp]-[i^-C)q-vp] 

7o Wo 'J {27rY [p - + ({1 - C)q^ + r]{l - r])p'^f' ^ ^ 

As a consequence of eq. (|37|), the term involving T^^^ will always carry an explicit factor of 

2.4 Total Energy 

We can obtain the total energy simply by integrating the energy density in eq. ([T6|) , 

E[a] = / e{r)dr. (44) 



Since both the t integral and the sum over channels in eq. (|1^) are absolutely convergent, 
we can interchange the order of integration and obtain 

E[a] = -Y.NJ -VW^^ dr[peitt,r)]^ + J2E^^l + EcT, (45) 
^ Jm n Jo -^^ 

where the total derivative term has integrated to zero. As explained in the previous subsec- 
tion, we have to subtract a sufficient number of Born approximations to the local spectral 
density pi{it,r) to render the t integral convergent. These subtractions are then added back 
in as the contribution to the total energy from the Feynman diagrams 



^=l ^=l 

which we then combine with the contribution from the counterterms 

EcT= / drecT{r) (47) 

JO 

to obtain a finite result. In practice, these terms can be more efficiently computed directly 
from the perturbation series of the total energy. For the first term in eq. (|45|) , we employ 
the relation 

2 dr[pe{k,r)]^ = — dr Ge{r,r,k)-Gf\r,r,k) =i—\nF,{k) (48) 

JO X Jo drv 



which we prove in Appendix A for all k in the upper half-plane, Im > 0. Using eq. (|48| ) on 
the imaginary axis k = it, we find in particular 

2 Hdr [p{zt,r)], = ^\nF,{zt) = Un hm -MM = Ung,{tt,^). (49) 
JO dt dt r^o Wi{ztr) dt 
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Finally, we introduce the real function f3g{t,r) = In gi{it,r) and write 

dt r-^ ^ d ^ 



E 



o- 



Y^N.j ' ^VW^j [/3,(t, 0)]^ + E + EcT 

1=1 



dt 



t 



N 



271 



CT 



i=l 



where Pi is determined by the differential equation 



(50) 



(51) 



with the boundary conditions linir^oo /5(t, ?") = lim^^oo /5'(^5 ''^) = 0. To compute [/?£(t,0)]^, 
we must subtract the first N Born terms, which we compute by iterating the differential 
equation (|5T| ) according to the expansion of Pe(t,r) in powers of the interaction a. 

The formulation (|50D is well suited for numerical evaluation and we will use it exclusively 
to compute the total energy. To make contact with previous work |T5|, however, it is in- 
structive to return to the real axis. We rotate the contour in eq. (|50D to obtain an integral 
over the whole real axis, picking up the discrete contributions from the bound state poles. 
When k approaches the real axis from above, we have 



■ d , . d ^ ,r^/,M d5Ak) 

t \nFi{k)=z—\n\Fe{k)\ + -^ 

drv drZ (XK 



(52) 



where 5i{k) is the phase shift. The first term on the right-hand side is odd in k and does 
not contribute to the integral. Thus we have 



N 



+ J2Eg + EcT 



(53) 



i=l 



which was used in [^]. To understand the physical motivation for eq. (|53D , we introduce 
the difference between the free and interacting density of states. 



2k. 



— Im / (Gi{r,r,k + ie) -Gf\r,r,k + ie))dr = pi{k) - pf\k) 

TT JO ^ ^ 

which from eq. (HSI) and eq. (E2|) is related to the phase shift by 



Tc dk 



Thus we can write the energy as 

r°° dk 



E[a] = Y.Ne 



f-^ra[M^)]. + -E 



N 

E^: 

i=l 



{i) 
FD 



CT 



(54) 



(55) 



(56) 



and see that it is simply the sum of the zero-point energies ^u, including both discrete 
contributions from the bound states and an integral over the continuum, weighted by the 
density of states. 
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We can also see this relation directly by integrating eq. (^) over r, with Born terms 
subtracted and added back in as diagrams, 



lO+ie TT Jo j 2 ^ 

(57) 

where the total derivative term has integrated to zero and the bound state wavefunctions have 
integrated to unity. We can then use eq. ( [A.14| ) to see that eq. ( ^71 ) is equivalent to eq. (|56[). 



While this relation is useful for building an intuitive understanding of the expressions for the 



total energy and the energy density, eq. ( [57| ) is of no practical use in numerical calculations; 
it only converges because eq. (0) has introduced the standard ie prescription to control the 
oscillations of the wavefunctions at spatial infinity. Any direct numerical calculation using 



eq. (|57D will be hopelessly obscured by these oscillations. 

3. Examples in One Space Dimension 

To illustrate our method, we apply it first to simple Casimir problems in one space 
dimension. There are only two channels (symmetric and antisymmetric) in this case. The 
boundary condition in the symmetric channel, ^|a;=o = 0, requires slight modifications of 
the formalism developed so far. We study two examples: a background field consisting of 
(5-functions, which allows us to make contact with treatments already in the literature p, ^ , 
and a GauBian background, to demonstrate the efficient use of our methods for numerical 
computation. 

3.1 Green's Function 

In one spatial dimension, the S-matrix for a symmetric potential is diagonalized in a 
basis of symmetric and antisymmetric wavefunctions. In the antisymmetric channel, we can 
proceed with the general techniques described in Section 2. From eqs. ( pSj ) and ( P7[ ) we 
compute g{it,x) and h-{it,x) with ^ = 1. Note that g{it,x) is the same in both channels 
because the channels only differ by the boundary conditions at the origin, which do not 
affect g. At equal spatial arguments we have 

^ / ■.^ h_{it,x)g{it,x) 

G_ X, X, tt) = — — . 58 

g{tt,0) 

However, the boundary conditions in the symmetric channel, 

0+(O) = 1 and 0'(O) = O, (59) 



are not in the class defined by eq. (pOl). Constructing the Jost solution f{k,x) with the 
boundary condition limx^oo f{k,x)e~^^^ = 1 and working out the Wronskian for 0+ and / 
yields the Green's function in the symmetric channel 

G+{x,y,k) = 'flkO) ' 
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where f'{k,0) = dxf{k,x)\x=o- By symmetry we have G^{—x,—y,k) = G+{x,y, k). We 
factor out the asymptotic behavior, 

(j)+{k,x) = e-'''''h+{k,x) (61) 

and obtain 

- hl{it, x) - 2th'^{it, x) + a{x)h+{it, x) = . (62) 

which is the same differential equation obeyed by h^{it,x). The boundary conditions in 
eq. (|59D then become 

h+(it,0) = l and h'^{it,0) = -t . (63) 

With this solution, we can now construct the Green's function in the symmetric channel. To 
compute densities it again suffices to consider x = y, where we have 

^ / ■.^ g{it,x)h+{it,x) 

G+(x,x,it) = — ; ^ 64 

and the complete Green's function is the sum G — G-\- ~h G ^ giving 

/■+ N .9{it,x)h_{it,x) g{it,x)h+{it,x) 

pUt,x)=t z ^ z ^ z r— . 65 

The Born series is then obtained by iteration of the differential equations (pS]), (p7|), and 

dH). 

3.2 Total Energy 

The contribution of the symmetric channel to the total energy cannot be cast in the 
general form of eq. ( pOD because, as for the Green's function, the boundary condition for the 
corresponding wavefunction is not of the type of eq. (^). Since the derivative of 0+ vanishes 
at X = 0, the Jost function becomes 

.,^ , f{k,x) , o'(fc,0) , , 

a;^0 ik ik 

which is normalized so that F^{k) = 1 when a = 0. Otherwise, the same analyticity 
arguments apply as in channels whose wavefunctions vanish at r = 0. In particular, the k 
integration contour can be closed around the discontinuity along the positive imaginary axis. 
We obtain 



E[a] 



dt t 



271 y/t 



2 _ ^2 



\ng{it,0) +\n[g{it,0) -jg'{zt,0) 



+ Ei'^ + EcT , (67) 



where as before the subscript indicates that the first Born term has been removed from 
the expression in brackets. One subtraction suffices in one dimension. By our no-tadpole 
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renormalization condition, the counterterm cancels the tadpole diagram identically, so that 

+ ^CT = 0. 

3.3 Delta— Function Background 

We first consider the case of two (5-functions separated by a distance 2a, 

as{x) = X [6{x + a) + 6{x - a)] . (68) 

At each 5-function, the derivative of the field jumps by For modes with /c -C A, the 
wavefunction il){k,x) vanishes at ±a, so that the limit A ^ cxd gives a Dirichlet boundary 
condition at ±a. Since the counterterms are local functionals of the background field, they 
are identically zero in regions where the background fields vanish. Hence in a renormalizable 
theory, local observables such as the energy density cannot have ultraviolet divergences in 
these regions, as we will see explicitly below. 

By integration of the respective differential equations we obtain 



g{it,x) 
h^(it, x) 



h-{it, x) = — 




\x\ > a 

< Ixl < a 



1 + e 



-2ta 



-2tx 



-2t{x-a) 



^-'ita _|_ g-2ta; _ ^-2t{x-a) 



\x\ > a 
< |x| < a 

|x| > a 

< Ixl < a 



(69) 



It is clear that g, and are all well-behaved for t > 0. We substitute g{it,0) and its 
derivative into eq. (p7|) and find the total energy 



^2 (a. A) 



dt 



1 



27rv^ 



tin 



1 + 



A A^ 

t ^ 4t2 



-Ata 



A 



(70) 



The last term in curly brackets is the = 1 Born subtraction. Carrying out a similar 
calculation for a single (5-function, we obtain 



^ ^ 2t 



27C 

A — mvr 



47r 



+ 



m 
2^ 



-;— 2 arccos 

4m*^ 2m 



Am? -'- ( 2m 



Am? 



A < 2m 
A > 2m . 



(71) 
(72) 



From eq. ([7^ ) we can verify that Ei and E2 obey the consistency conditions lim £"2 (a. A) 
2Ei(A) and limE2(a, A) = £i(2A). 

a— >0 
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We observe that in the hmit of large A the total energy associated with a single (5-function 
approaches minus infinity as —A In A, which cannot be canceled by any available counterterm. 
Thus in this limit the energy is infinite. In this sense the Dirichlet-Casimir problem is not 
well defined in the context of renormalizahle quantum field theory. However, we can compute 
the force between two 5-functions for finite A and subsequently take the A ^ oo limit, 



\ I , . n/r,„\ / _ /To ^ilMa 1^ ^ > 



A-»oo (9(2a) Jm -K \ft 



(e 



which is equal to the result found using boundary conditions ||^. The massless limit does not 
exist for the vacuum polarization energy in eq. ([7^), even for finite A. This is to be expected 
since massless scalar field theories are infrared divergent in 1 + 1 dimensions. However, it 
exists for the force, yielding the well-known result 

limi^(a) = --^. (74) 

Next we turn to the energy density. In this Section, we consider the energy density 
e(x) = (Too) without factoring out the surface area factor as we did in the previous Section 
(which would just give a factor of 2 here). The resulting energy density e{x) is normalized 
such that /f^ e{x)dx gives the total energy. It is convenient to decompose the energy density 

into e(x), the integral over the subtracted local spectral density, and ep^^ix) + ecT{x), the 
contribution of the Feynman diagram plus counterterms. 

First we study e(a;). From the radial functions in eq. (^), we find the local spectral 
density ClS db function of x and k = it. 



A 



2t - A + (2t + A)e 



4ia 



,2t{a-\x\) 



A2 - {2t + A)2e4*« 



(2t + A)e^*"cosh(2ta;) - A 



< a; < a , 



where only the free local spectral density has been subtracted. For 7^ a the local spectral 
density decays exponentially as t increases, so it yields a finite energy density without any 
subtractions, as we expected since the counterterm vanishes away from x = ±a. 

To prepare for more general problems involving smooth potentials, where the energy will 
diverge at all x, we carry out the full Born subtraction procedure even though in this case the 
Born subtraction and the diagram contribution are both finite and equal for all 7^ a, and 
so their contributions cancel. We will thus obtain the energy density everywhere including 
the points x = ±a. The first Born approximation is given by 



p(^\it,x) = -^{ (76) 



2t 



and we obtain 



-/ u I dt r]{it,x) 

e{x,a,X) = / (77) 
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where 

{2t + A)2e4*'^ - A' I ^(^2 _ ^2) ^ ^(2t + A)e2*« - Ae'^*"] cosh(2te) < |x| < a . 



-FD 



For the Feynman diagram and the counterterm, we use eq. (^if) and spht off a (5-function 
contribution, 

e[.\^(x, a, A) + ecT(a;, a, A) = -^5{\x\ - a) 

A f°° TTL^ 

dq cos{qa) cos{qx) / dC^—— — ^ . (79) 

Tc^ Jo Jo + C(l — OQ' 

We recast the second term in a form which makes its identification with the first Born 
contribution, eq. ([7BD, manifest, 

, , r/i I \ 4m^A f°° , cos(ga) cos(gx) . , _i g 
(x,a. A) + ecT(a;,a, A) = — - a) — / rfg . . smh — 

which confirms that our Feynman diagram representation of the 0{X) contribution to 
e(x, a. A) agrees with the k integral representation for values of x where the counterterms 
vanish. 

In Fig. ^we display the energy density for two delta functions separated by 2a, 62 (x, a, A) 
as functions of x. Note that there are 5-function contributions to €2(0;, a. A) at a; = ±a not 
drawn in the figure. For reference we also show the sum of the energy densities for delta 
functions at x = a and x = —a separately, which we denote ei(x =F a. A). The difference 
between e2(x, a. A) and the sum ei{x — a. A) + ei(x + a. A) will result in a force between the 
two points. 

As A — > cxD, we expect to obtain the same energy density as if we had imposed the 
boundary conditions 0(— a) = 0(a) = from the outset. No subtraction is necessary away 
from X = ±a, so we can study the energy density in this region directly from [p]q, given in 
eq. (P, 

1 f°° dt — m? + rn^e^^"- cosh(2tx) 
hm e2 X = / 2 — for \x\ < a . 81 

To compare with Ref. we separate the x-dependent part and transform the t integral 
into a sum over the poles of the function e^*'^/(e'^'^* — 1) = l/(2sinh2ta) at t = m7r/2a. 



I dt t^-m^ m m2 „ cos ^(|x 

lim eoix) = / =— — > — , 

Try™ ./W^^e^-^-l 8a 4a . /7^^ ^2 



for Ixl < a . (82) 
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Figure 3: Left panel: The total energy density e2(x,a, A) and the contributions e{x,a,X) and 
epl){x, a, A), defined in eqs. ([77|) and (|80|). Note that the displayed energy densities are supplemented 
by (5-function contributions located at x = a and x = —a. Right panel: The total density e2(x, o, A) 
compared to the sum €i{x — a, A) + ei(x + a, A) of the energy densities of two isolated (5-functions. 
The parameters are A = 3m and a = 0.2/m. 



This expression agrees with the Dirichlet result given in Ref. 0. 

This simple example allows us to study the nature of the divergence in the Dirichlet 
Casimir energy in the limit A — > oo. For any fixed x away from the 5-functions, the energy 
density reaches a finite limit as X oo, but the limit is approached nonuniformly. This 
can be seen clearly in Fig. ^, where we display the energy density for large couplings and 
compare it to the limiting function, eq. (0). The closer x approaches to a, the slower the 
energy density converges to its limiting form. Thus the energy density away from the source 
on the boundary remains finite and well defined in the boundary condition limit even though 
the total energy diverges. 

3.4 Gaufiian Background Field 

In this subsection we demonstrate how to calculate the energy density associated with 
a GauBian background using our method. This problem is a warmup for the next Section, 
where we will use a tall, narrow Gaufiian background as an approximation to the 5-function in 
higher dimensions, where the energy of a delta function diverges even at finite coupling. We 
choose a Gaufiian background that approximates the 5-function potential from the previous 
subsection. 



A r (x-a)^ (x + af 



(83) 



in the limit w ^ where A is chosen so that the area under the Gaufiian is fixed to 2A. We 
substitute this background field into the differential equations (^Sf), (|27| ) and (^) and obtain 
the wavefunctions g, /i_ and numerically. As promised, these functions are well-behaved 
on the positive imaginary k axis. At large momenta t, we expect /i_ {it, x) to approach the 
free solution, which decays asymptotically as (1 + e~*'^'')/2t. This behavior can be seen 
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X/m 

Figure 4: Ratio of the energy density for two (5-functions with finite A to the energy density for 
A ^ cxD with a = 0.2/m, plotted as a function of A for several values of A = (a — a;)/a. 




Figure 5: Plots of g{it,x)/g{it,0) and h-{it,x) as functions of the imaginary momentum t at 
several positions x. The parameters are a = 0.3/m, w = 0.03/m and A = 3m. 



21 



Figure 6: The Born subtracted weight function r]{it,x) for the Gaui3ian background field, eq. (|8^), 
as function of the imaginary momentum t at several positions x. The parameters are a = 0.3/m, 
w = 0.03/m- and A = 3m. 



clearly in Fig|^. By the same argument, we expect the ratio g{it,x)/g{it,0) to approach one 
at large t, which is confirmed in Fig ^. The functions required for the symmetric channel 
behave similarly. The Born subtracted continuum contribution to the energy density can be 
written in the form of eq. (|77D , where now ri{it, x) is computed numerically from g{it, x) and 
h±{it, x). In Fig. ^ we plot ri{it, x) for the GauBian background for several values of x. Since 
a is nonzero everywhere, t] falls like at large t for all x. We see that the contribution 
for I X lb a I ^ If , where the potential a is very small, is much smaller than the contribution 
from X ~ a, where the potential is peaked. 

The energy density we obtain for a tall, narrow Gaufiian is shown in Fig. together with 
the result from the previous subsection for the ^-function background. The latter case also 
has a (5-function contribution to the energy, which is not shown. For |x — a| > 3w the energy 
densities for these two background fields are indistinguishable. For |x — a| < 3w the energy 
density for the GauBian background becomes dominated by a strongly peaked contribution 
proportional to cr(x) (approximating the (5-function contribution). 

Finally, as a check on our numerical calculation, we have verified that the total energy 
computed directly from eq. ( p7| ) agrees with the spatial integral of the energy density. 

4. An Example in Two Space Dimensions 

Our numerical methods become particularly valuable as we move to higher dimensions. 
Our renormalization procedure and numerical methods allow us to study a sequence of 
sources and examine the behavior of the energy and energy density as one attempts to take 
the boundary condition limit. In two dimensions we still only need the renormalization 
counterterm proportional to a{x). To go on to three dimensions would require an additional 
Born subtraction and a counterterm proportional to 



We begin by discussing some peculiarities of the differential equations (^5]) and ( pTf ) for the 
auxiliary radial functions ge and hi necessary to handle the s-wave in two dimensions. Then 
we describe the computation of the total energy and energy density for a GauBian background 
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Figure 7: Comparison of the energy densities for a narrow Gaui3ian, ac with w = 0.003/m, and 
a (5-function background, as, located at a = 0.3/m. The couphng is A = 3m. Note that we have 
omitted the 5{x — a) piece for ag. Its coefficient is such that when it is added to the dotted hue 
the areas under both curves are equah 



of width w centered on a circle of radius a, which illustrates our general conclusions about 
the nature of the Casimir problem. 

4.1 Boundary Conditions for Radial P\inctions 

In two space dimensions the angular momentum i can be restricted to ^ > with degen- 
eracy Ni = 2 for i ^ and Ni = 1 for £ = 0.0 In each channel with £ > 0, we can solve the 
radial problem by following the steps described in Section 2.2. Specifically, u = £ — 1 + ^ = i 
and 

e.W = -#ln(v^i^.(^)) = ^ + %# ii>0), (84) 



dz ^ ' z Ke(z) 

where K^i = Ki. These functions are smooth for positive arguments z > and quickly 
approach unity for z ^ 1. For i > they have a simple pole at 2; = 0, C,e{z) — {£ — \)/z. 
However, the functions and are regular near r = 0, 

g,{it,r) = g,{it,Q) + 0{r^) h,{tt,r) = r + 0{r^) (85) 



canceling the potential singularities at r = in eqs. (25) and (27) coming from ^e{z). 

For £ = it appears from eqs. ( p3l) and ( p4|) that both ho and the Green's function Go 
become ill-defined for u = i = 0. Ultimately, the source of the problem can be traced to the 
subleading logarithmic term in the Bessel functions, 

woikr) = i^l^Hl,'\kr) and ^o(^) = - + O(^) • (86) 



^In the formula for Nt, given after eq. we first have to take £ — and subsequently n — > 2 for the 
proper treatment of the singularity in r(n + £ — 2) for this case p3|. 
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The subleading term shows up as a singularity in Coi^)^ while it is harmless in Co{z) itself. 
The solution is to integrate the differential equation (^) for qq from r = oo down to a point 
arbitrarily close to r = and then get go{it, 0) from eq. (pS]). 

For the regular solution, we simply omit the normalization l/2t/ in the ansatz (|23|) and 

set 

Vkhoik,r) 

Mf^^ r) = TT-^ ■ i^V 



Then eq. (^) becomes 



Go{r, r,k) 



Wo{kr) 
go{k,r)ho{k,r) 



9o{k,0) ■ ' ' 

The resulting differential equation for ho is still eq. (|27|) , but the boundary conditions at 
r = must be modified: It is straightforward to verify that 



hQ{it,r) = —rlntr and /iQ(zt, r) = — [1 + Intr] 



as 







^9) 



removesQ all small r singularities in eq. (P?]). We use eq. ( PD| ) to integrate the differential 
equation (p7|) numerically starting at an arbitrarily small but nonzero value r. 

It is then also straightforward, although tedious, to expand and hi in powers of the 
background field a according to eq. (^). We thus obtain the Born subtracted local spectral 
density, eq. (14), for each angular momentum quantum number i. 

The final computational ingredient in eq. (|T6|) is the renormalized contribution from the 
Feynman diagrams. In two space dimensions one Born subtraction suffices. Then we have to 
add back in the two Feynman diagrams in Fig. |^. In the no-tadpole renormalization scheme 
the local diagram, eq. (pHD, is exactly canceled by the counterterm and we only need to add 
back the two-point contribution, eq. (|1T|). The corresponding Feynman parameter integral 
is readily computed to be 



167r JO 



dpJo{pr)a{p) 



p + 



p 



— 4m p 

arctan 

2m 2m 



where 



a{p) = 27r / dr rjQ(pr)a(r). 



(90) 



(91) 



4.2 Numerical Results for a Gaufiian Background Field 

We now apply this method and discuss the numerical results for the total energy and 
energy density of a Gauf^ian background field in two dimensions. We consider 

cr(r) = AAe"^^ . (92) 

^The logarithmic behavior is actually expected from the free solution: t/q (tr)Kci (tr) — > In 2 — 7 — r ln(ir ) 
as r ^ 0, cf. eq. (|3l|). 
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Figure 8: Contributions Ei{X,a,w) to the total energy from various angular momentum channels 
£. The logarithmic plot on the right is consistent with the asymptotic behavior Eg ~ The 
parameters are a = 1.0/m, w = 0.1/m and A = 3m. 



This profile describes a ring of radius a and thickness w. The normalization A is chosen such 
that a{r)dr = A, giving cr(r) — > X6{r — a) as w 0. Sending the coupling strength A to 
infinity subsequently imposes the Dirichlet boundary condition on a circle of radius a. 

We first discuss the total energy E[a] = E{X,a,w) computed from eq. (|50D. Note that 
no Feynman diagram contribution has to be added to the total energy because in the no- 
tadpole scheme Ep^ + Eq^ = 0. Fig. ^ shows the total energy contribution Ei{\,a,w) 
of the angular momentum channel i, for a Gaufiian ring of radius a = 1.0/m and width 
w = 0.1/m. Note that the i = channel is suppressed relative to £ > by the degeneracy 
factor No/Ni = \. Unlike the case in three dimensions, however, higher angular momentum 
channels are not favored by a further increasing degeneracy, so the contributions decay 
rapidly with increasing i > 1. The largest contribution always comes from i = 1. The 
logarithmic plot in Fig. § confirms that the asymptotic region is reached quickly and that 
the decay of the individual contributions is consistent with Eg oc l/i"^ for i > 5. While this 
clearly shows the finiteness of the renormalized quantum energy (for fixed a and w), the 
convergence is too slow to allow us to sum the ^-series directly. Instead, we use a speedup 
technique related to Richardson's approach ||2^, which is particularly powerful in the present 



case, where the asymptotic behavior sets in quickly. We have checked that the results of this 
method are very robust and independent of the details of the extrapolation technique used. A 
sufficiently high precision in the energy and density calculation can then be achieved by using 
only approximately 30 channels. For details on this summation technique, see Appendix C. 

The computation of the total energy can be sped up even further by subtracting the 
second Born approximation in eq. (0) and adding back the corresponding second-order 
Feynman diagram. With this (over-) subtraction, the large £ behavior of E^{\,a,w) changes 
to The relevant Feynman diagram is easily computed, 

^fdM = ^^^ dpa{pf w:ciiin^ (93) 
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Figure 9: Contributions Ei to the total energy from various angular momentum channels i. The 
GauBian background is much narrower than in Fig. ^and the asymptotic decay Ei l/£^ is reached 
for much larger channels £ only. The parameters are a = 1.0/m, w = 0.008/m and A = 3m. 



where a denotes the Fourier transform, eq. (j^^. We have verified that this method gives 
the same result as before. 

As the background field approaches a 5-function we find by explicit numerical calculation 
that the total Casimir energy diverges. To illustrate the problem we consider a very narrow 
configuration with w = 0.008/m. As can be seen by comparing Figs. ^ and ^ Ei decreases 
more slowly with i when w is smaller. For any finite width w the energy contributions Ei 
eventually decrease steeply enough to yield a convergent sum and a finite total energy, but 
this asymptotic regime is only reached for much larger i when the width is small. This 
is a typical non-uniform behavior: For any fixed width, there is always a channel £o large 
enough such that the energies Ei decay as for £ > £q. On the other hand, for every 
fixed channel £ we can always find a width w small enough such that we are far from the 
asymptotic region. Eventually, as w ^ the asymptotic region is never reached and the 
total energy diverges. 

One can also understand this divergence in terms of the second Born approximation, 
which we separated out in the oversubtraction procedure above. For the singular background, 
a{r) = X6{r — a), this diagram contributes 



Epl\X6{r — a)] = / dpj^ (ap) aictan (94) 

8 JO 2m 



which diverges logarithmically at large p just as the sum over channels did without oversub- 
traction. Oversubtraction shifts the divergence from the £ sum into the momentum integral 
in the Feynman diagrams. Whichever choice we make, the divergence cannot be renormal- 
ized away, despite the claims of |T^; it is a physical effect. Its implications for the usual 
Casimir calculations will be studied elsewhere [0. 

Next we examine the radial energy density e(r). As a numerical check, we have verified 



that its integral over r agrees with the total energy computed directly from eq. (0). In Fig. |T0 
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Figure 10: Left panel: Contributions of various angular momentum channels to the Born sub- 
tracted energy density e(r). Note that we have multiplied the i = contribution by a factor of 
2 to better display the secular decrease with i. Right panel: The Feynman diagram contribution 
eFD{r). The parameters for both panels are a = 1.0/m, w = 0.1/m and A = 3m. See also Fig. ^ 
for the total energy density e(r) found by combining the components displayed here. 



we display the contributions of various angular momentum channels to the Born subtracted 
energy density, e{r), for a relatively wide profile, w = 0.1/a. Generically the contributions 
to the energy density from various angular momentum channels vary in magnitude but 
have similar radial shapes. As for the total energy, the i = channel is suppressed by a 
degeneracy factor 1/2; the largest contribution comes from i = 1 with all higher channels 
decaying quickly with increasing i. For narrower widths the decay with i is slower, as we 
expect from the preceding discussion of the total energy. The smooth asymptotic behavior 
allows us again to use an extrapolation method in order to speed up the sum over channels. 

Figure |10| also shows the renormalized Feynman diagram contribution in the right panel. 
Since the integral over all radii r is zero, the peak at the location of the background profile 
must be compensated by negative densities on both sides. The sign and width of the peak 
essentially reproduce the GauBian background itself, while the similar peak in the density 
e(r) has the opposite sign and reduces the overall magnitude of the density. As can be seen 
from the total energy density e(r) in Fig. |ll], combining all the contributions still yields a 
pronounced peak at the location of the background. As we decrease the width of <j{r), the 
peak increases in magnitude and becomes narrower. In the (5-function limit it develops a 
non-integrable singularity, unlike the case in one dimension. 

Figs. |Tl| and |1^ display the dependence of the energy density for various values of r on 
the width of the GauBian background field. We see the non-uniform approach to the limit 
w ^ 0. In particular, it is clear that the energy density approaches a finite limit at any fixed 
r a.s w —>■ 0, even though the limiting function diverges as r — > a. Sufficiently far away from 
the location of the GauBian, around |r — a| > 3w, the energy density quickly approaches 
the limiting form, corresponding to that of a 5-function background. In contrast, at r = a, 
no such limiting value exits. This behavior can also be seen in Fig. where we display 
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Figure 11: The energy density in units of m? for the GauBian background field located at a = 0.1/m 
with strength A = 3.0m for various values of the width w. The left panel shows the energy density 
for all radii. The right panel focuses on the region where the densities can be seen to converge 
toward the limiting form. 



the energy density at various values of r as function of the width. The limiting form is 
approached more rapidly the farther one is away from the GauBian peak. 

The Dirichlet boundary condition is achieved only by taking a{r) XS{r — a), followed by 
the limit A — >■ cx). However, since the renormalized Casimir energy diverges as a approaches 
a (5-function, it is impossible to define the Casimir energy for any sharp background field, 
much less one with infinite coupling. Thus the Dirichlet-Casimir energy of a circle in two 
dimensions does not exist. But all is not lost: the energy density away from the source 
approaches a finite limit as the background becomes sharp. 

5. Summary and Outlook 

In this paper we have developed an efficient method to compute one-loop energies and 
energy densities in renormalizable quantum field theories. Starting from a Green's function 
formalism, we employ Born subtractions to implement definite renormalization conditions, 
and express the result in a form that allows for direct, efficient calculation after analytic 
continuation to the imaginary k axis. This method allows us to study the limit in which 
smooth background fields approach idealized boundary conditions, illuminating subtleties of 
the process by which standard Casimir calculations are realized in renormalized quantum 
field theory. We find that although the Casimir force between rigid bodies and energy 
densities away from the boundaries are well-defined observables, quantities like the Casimir 
energy have divergences that cannot consistently be renormalized. These divergences have 
physical consequences in examples such as the stress on a circle, where we have to compare 
two configurations for which the divergent contributions differ. In these cases, the physical 
effects depend strongly on the details of the source. 

We have also shown rigorously that the total energy, as obtained from the integrated 
energy density, can be expressed as the integral over single particle energies weighted by the 
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Figure 12: The energy density of the GauBian background field located at a = 1.0/m and with 
strength A = 3.0m. The figure shows the energy densities at various values of r as functions of the 
width, w. The densities are scaled to e(r) at ti; = 0.1/m. 



change in the density of states, which is given by the derivative of the phase shifts. 

There are many possible apphcations of our method. It can be generahzed to fermions 
and to situations where the energy spectrum is asymmetric about zero. In the latter case we 
would have to add the different contributions from the two sides of the cut in Fig. |l]. In such 
systems the background field can carry nonzero charge density, since the contribution to the 
charge density from positive and negative energy modes will no longer cancel. It could also 
be calculated with our approach. We could consider the electromagnetic field in the presence 
of conductors by decomposing it into a scalar field obeying Neumann boundary conditions 
plus a scalar field obeying Dirichlet boundary conditions. We have seen in this paper how 
to implement Dirichlet boundary conditions with a background potential, and the Neumann 



boundary condition can be treated in a similar way [^. Therefore we are equipped to 



reconsider problems such as the surface tension of a conducting sphere in QED [Q. Our 
calculational method is well suited to soliton calculations, generahzing beyond the case 
of refiectionless potentials. Finally, in cases like the Gross-Neveu or Nambu-Jona-Lasinio 
models, where one introduces the background potential as a non-dynamical auxiliary field 
coupled to a dynamical field, one can search for solitons as self-consistent solutions to the 



auxiliary field equations of motion ||2^, |29[ . These equations of motion are obtained by first 
integrating out the dynamical field, leaving an effective action for the auxiliary field. This 
procedure generates expressions involving the full Green's function similar to those we have 
considered here, which could be computed efficiently with our methods. 
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Appendix A: Density of States and the Jost Function 

\) between the integral over r of the 



In this Appendix, we demonstrate the relation 
local spectral density and the Jost function 



2k 

dr [pi{k,r)]Q = — 

I JO 



dr 



Ge{r, r, k) - Cf (r, r, k)\ = In F,{k) 



d 



(A.l) 



which is valid everywhere in the upper half-plane \mk > 0. This formula can be used both 
to establish an efficient computation of the total energy on the imaginary axis, and also to 
relate the density of states for real k to the phase shift. 

We start by differentiating the Wronskian of the Jost solution, /^(fc, r), and the regular 
solution, (f)i{k',r) [pO|] , 



[f,{k, r) , Uk\ r)] = - k") fe{k, r)Mk', r) 



(A.2) 



where all quantities in this relation are analytic for Imfc > 0. We integrate both sides from 
r = Q to r = R. Since the regular solution (l)£{k,r) becomes fc-independent at small r, we 
can compute the boundary term at r = by replacing k' with k and using the standard 
Wronskian, 

(A.3) 



W[Mk,r),Mk,r)] = {-k)-^~''F,{k) 



givmg 



W[feik,R),Mk',R)] = i-k)-^-'' Feik) 



J2 



R 



dr ft,{k,r)(f)i{k',r) 



Next we differentiate with respect to k, set k' = k, and use the representation 
Green's function to obtain 



:-kr 

Fi{k) 



-w 



Mk,R),Mk,R) 



Feik) 



R 



2 +^^ + 2k I drGiir,r,k), 







(A.4) 
for the 

(A.5) 



where fi{k,R) = -^fi{k,R) and Ff(fc) = -^Fi(k). To eliminate the first term on the right- 
hand side, we subtract the same equation for the non-interacting case, giving 



:-ky 



Fiik) 



-w 



Mk,R),Mk,R)] - {-kr-'^w[f}'\k,R),<pP{k,R) 

Gi{r, r, k) — Gf^ (r, r, k) 



^ ' + 2k dr 



Fe{k) 



(A.6) 



To complete the proof, we have to show that the left hand side of eq. (|A.6| ) vanishes as 
R oo. To see this, we write the boundary condition (|18D for the Jost solution in the form 



Mk,R)=w,{kR) \l + 0{R-'] 



R oo 



(A.7) 
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which can also be inferred from the integral equation obeyed by fe{k, r). Differentiating with 
respect to k and using the asymptotics of the free Jost solution Wi{kK), 



Wi{kR) 



d_ 

dk 



Wi{kR) = iR we{kR) 1 + 0{R- 



it is easy to show the asymptotic behavior 

Mk, R) = iR feik, R) [l + OiR' 
The first term on the left-hand side of eq. ( [A .61 ) can thus be estimated by 

i-ky ' 



(A.8) 



Feik) 



w[Mk,R),Mk,R)] 



iR 



-ky 



-W[feik,R),Mk,R)]-t 



'-ky 



Feik) 

t[R + Ge{R,R,k)] \l + 0{R-^) 



Feik) 



feik,R)Mk,r)} 1 + 0{R 



(A.9) 



where we have used the Wronskian of fi and (pi and the definition of the Green's function, 
eq. (^). Subtracting the analogous equation in the free case, the term proportional to R 
drops out and we are left with 



FfXk) 



-w 



Mk,R),Mk,R)] - {-kr-'^w[f}'\k,R),<f)P{k,R) 
Gi{R,R,k)-Gf\R,R,k)] \l + 0{R-^ 



(A. 10) 



We estimate the large-i? behavior of the difference ^i{k,R) = Gi{R, R,k) — Gf\R,R,k) 
from eqs. (|A.9|) and (|A.6|) , 

-tAe{k,R)[l + 0{R-')] = -^ + 2k drAe{k,r). (A.ll) 

Fe[k) Jo 

It is then straightforward to solve for A£{k,R) ^ i[Gg{k) + 0{1)R\ exp{2ikR) at large 
R, where Gi{k) is an /^-independent integration constant. The left-hand side vanishes 
exponentially for all ImA; > as i? ^ oo, so in this limit we obtain eq. (|A.1|) , which 
completes the proof. 

It is also instructive to explore the consequences of this result on the real axis. If we let k 
approach the real axis from above, we see that to leading order in R~^, the integral equation 
(|A.11| ) allows for a solution going like Rexp[2ikR], which violates the bounds in eq. (|30|) if 
k is real. Thus the coefficient of this term must actually be zero, leaving 



Mk) 
Fiik) 



+ 2k I dr Gi{r,r,k) -Gf\r,r,k) = Gi{k) exp{2ikR) 1 + 0{R-^) (A.12 
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which oscillates as i? — * oo when k is real. As is typical for continuum problems, we must 
specify that the limit where k becomes real is taken after the integral over r to eliminate the 
contribution from these oscillations at the upper limit of integration. 

On the real axis, we can relate the Jost function to the phase shift by eq. (^2|). Taking 
the imaginary part of eq. (|A.1| ) and using eq. (Q) yields the relationship between the density 
of states and the phase shift, 

= ^Im r (Ge{r, r,k + ze) - (r, r,k + ze)) dr = p,{k) - pf\k) . (A.13) 

TT ak TT Jo ^ ' 

We can also rewrite eq. ( |A.13| ) as 

- rrfrte(A;,r)^,(A;,r)-#*(A;,r)#(fc,r)) (A.14) 

where the momentum on the left-hand side is understood to be defined with the ie pre- 
scription to eliminate the contribution from the oscillations of the wavefunctions at spatial 
infinity. 

Appendix B: Two Jost functions in One Dimension 

In this Appendix, we discuss a technique special to one dimension, which applies even 
when the background field is not symmetric. Scattering in one dimension can be described 
in terms of two Jost solutions f±{k,x), which obey the boundary conditions 



X— >±oo 

The Green's function then is 



lim f±{k,x)e^'^'' = 1. (B.r 

X— >±oo 



Gix,y,k) = ^^/+(fc,x>)/_(A;,x<) , (B.2) 

where T{k) is the transmission coefficient. For numerical computations it is convenient to 
factor exponentials out of f±{k,x). We define 

9±ik,x) = e^^'^f±ik,x), (B.3) 

which obey the differential equations 

- g'^{it, x) ± 2tg'^{it, x) + a{x)g±{it, x) = (B.4) 

with the boundary conditions 

lim g±(it,x) = l and lim g':(it,x) = 0, (B.5) 

and are bounded in the upper half /c-plane. The Green's function for equal arguments then 
reads 

G{x, X, it) = -^a+i^t, x)g-{it, x) . (B.6) 
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The transmission coefficient is most conveniently computed using the Wronskian for 
yielding 

G(x X it) = 9+{it,x)g-{it,x) 

' ^ g'+{tt, 0)g4tt, 0) - g'_{tt, 0)g+{tt, 0) + 2tg+{tt, 0)g.{tt, 0) ' ^ ' ' 

for the Green's function at equal spatial arguments. In this form, we can compute the Green's 
function without encountering any oscillatory or exponentially growing contributions. 

Appendix C: Slowly Convergent Series 

In the examples studied in the main text, the sum over angular momentum channels i is 
slowly converging, typically going as l/i"^ for large i. From a numerical point of view, such 
a decay is too slow for a direct summation as it requires to compute 0(1000) channels to 
get a relative precision of 1%. Even worse, the narrow GauBian proffies relevant for Casimir 
problems have a tendency to fall like 1/i before they eventually turn over to a behavior 
at channels with large angular momenta, £ > £o ^ 1- 

To sum such series efficiently using only a small number of channels, we employ a tech- 

oo 

nique related to Richardson's method p3. Given a slowly converging sum J2 CLe, "we define 
a function f(t) at the particular points tn by 

n 2 

f{Q = J2ae, ^« = ^T (n = 0,l,2,....) (C.l) 

which gives the value of the full infinite series in the limit t — >■ 0. We assume that /(O) can be 
approximated by a smooth extrapolation using the points t„, which accumulate around t = 0. 
Thus given values Oq, . . . , oat from the ffist + 1 channels, we construct an interpolating 
function / of a fixed order z/ through the last u + 1 points t^jt^^i, . . . ,t]\f_i, (for the ffist 
few channels N < u, so in those cases we take u = N). The extrapolation to t = gives an 
approximation /(O) for the sum of the series and an error estimate based on the assumption 
of a smooth underlying function f{t). While new channels are computed, we monitor /(O) 
and terminate if we have a satisfying error estimate. For our purposes, the best results are 
obtained from a rational function (or Fade) approximation, since the resulting f{t) has the 
tendency to anticipate the smooth behavior of the underlying f{t). To illustrate the method, 
let us look at the example 

which imitates the behavior of our typical sums over channels: It starts with ae ~ 1/i and 
turns into the converging l/i"^ behavior for i > i^. In the table below, we show the results of 
our method with the standard parameters used in the main text. We also give the number 
of channels actually summed and compare with the results obtained by direct summation of 
the same number of channels. 

As we can see from the table, our method improves the accuracy dramatically. For small 
turning points io, the ffist 10 channels are sufficient to find the exact value. As we go to 
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Table C.l: Results of our extrapolation method as compared to direct summation for the test 
series eq. ( p.2D . 



4 


5 


10 


50 


500 


1000 


5000 


10000 


exact 


2.2833 


2.9290 


4.4992 


6.7928 


7.4855 


9.0945 


9.7875 


Interpol. 


2.2833 


2.9290 


4.4949 


6.7385 


7.3077 


8.0560 


8.2627 


error [%] 


2.2 ■ 10-*^ 


2.38 ■ 10-5 


0.095 


0.80 


2.37 


11.4 


15.6 


direct 


1.8941 


2.2602 


2.7483 


3.9672 


4.1653 


4.2706 


4.3684 


error [%] 


17.0 


22.8 


38.9 


41.6 


44.4 


53.0 


55.3 


channels 


10 


10 


10 


31 


37 


40 


44 



higher it becomes more difficult to predict the exact result from the lowest terms which 
do not yet see the asymptotic decay. Even then, with ^Q as large as 500, we can still achieve 
a precision of < 1% using the ffist (9(30) channels only. The direct summation is off by 
more than 40% in this case. Eventually, for very large £0 — 10000, the error of the Pade 
approximation increases up to 15%. This is not too surprising since there is no trace of the 
eventual l/£^ decay in the ffist 44 channels summed. In order to predict the correct result 
within 1% for the case £0 = 10000, our method would in fact have to use (9(400) channels 
(but the direct summation would still be off by > 30% in this case). 
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